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Abstract 

The problem of vibration suppression in helicopter rotors is the 
main topic of this paper. We revisit a strategy for the stabilization 
of such systems based on delayed feedback. By means of an improved 
analysis, made possible by recently developed mathematical tools, we 
present insights on how the delay can be used as stabilization tool for 
this problem. The results are supported with numerical simulations. 
For simplicity, only the case of hovering flight is considered. 


1 Introduction 


The main rotor is the principal source of vibrations in a helicopter and the 
reduction of these vibrations is important to increase the life of the com¬ 
ponents of the airframe as well as the comfort of passengers and crew. A 
considerable research effort has been devoted to the study of active control 
systems for helicopter rotors, in order to guarantee the stability of its mo¬ 
tion. Most works on this topic fall in one of two categories: Higher Harmonic 
Control (HHC) or Individual Blade Control (IBC). Works in the HHC cate¬ 


gory, like those of Hall and Wereley (1992); Mura and Lovera (2014); Mura 


et al. (2014) consider the control forces to be applied to the swashplate by 


means of actuators. The IBC works (Bittanti and Cuzzola 

2002 

Arcara 

et al. 2000, 

King et al. 

2014 

Shen et al. 

2006 

) consider that each blade can 


be actuated independently. The main difference between the two approaches 
is the number of available degrees of freedom for the control action. Many 
control algorithms, nevertheless, can be implemented with either actuation 


technique. Friedmann and Millott (1995) compare several different control 


methods in both categories. A more recent survey of the developments in 


these topics was presented by Friedmann (2014). 
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The present paper revisits a control strategy presented by [Krodkiewski 
and Faragher (2000), which can be implemented in both HHC and IBC 


frameworks. This control logic uses delayed feedback to improve the stability 
of the motion of the helicopter rotor. Delayed feedback has been previously 


used for the stabilization of periodic motion by Pyragas (1992) and for the 


suppression of vibrations by Olgac and Holm-Hansen (1994). Using novel 


analysis techniques, we challenge the published results and show a corrected 
stability region in the parametric space. 


Although the original works of Faragher (1996) and Krodkiewski and 


Faragher (2000) consider both the cases of hovering and horizontal flight, we 


restrict the scope of this work to the case of vertical motion only. In this way, 
we can apply tools for the stability analysis of Linear Time Invariant (LTI) 
systems affected by delays. The case of horizontal flight, and the study of 
the Linear Time Periodic (LTP) systems with delays, is a matter of further 
research. 

The paper is organized as follows. Section [2] introduces the dynamic 
model of the rotor in hovering flight, borrowed from Bramwell et al. (2001). 


Section [3] presents the instabilities that may appear in the operation of an 
uncontrolled rotor. Section [4] introduces the control law and the stability 
analysis. Some concluding remarks are finally presented in section [5] 


2 Dynamics of the Rotor 

A helicopter rotor blade is mounted on a set of hinges which allow three 
angular degrees of freedom. A typical arrangement of the hinges is shown 
in figure [l] The flapping motion is defined as an up and down rotation in 
a plane which contains both the blade and the shaft. The flapping angle, 
j3, is defined as positive when the blade moves upwards. A flapping blade 
rotating at high speeds is subject to large Coriolis moments in the plane of 
rotation, and the lag hinge is introduced to alleviate these moments. This 
hinge allows the motion of the blade in the same plane of rotation. The 
lagging angle, £, is positive when the blade moves in the same sense of the 
shaft rotation. The pitching or feathering motion, denoted by 9, is a rotation 
about an axis parallel to the blade span. 

The pitch-flap flutter is an unstable behavior produced by the coupling of 
the pitching and flapping motions of the rotor. It can lead to very high oscil¬ 
latory loads in the pitch control mechanism. The following dynamic model, 


adapted by Krodkiewski and Faragher (2000) from the work of Bramwell 


et al. (2001), describes the coupling between 6 and ft for an uncontrolled 
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rotor in hovering flight 
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with the inertia, damping, and stiffness matrices: 
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As it is customary in helicopter dynamics, the independent variable in ([!]) is 
not the time t, but the azimuth angle ip = fit. Here fi is the angular speed 
of the rotor shaft, which is assumed constant during the operation of the 
aircraft. In (JTJ a dot over a variable denotes its derivative with respect to ip. 

It is important to state that in the derivation of ([!]) it is assumed that 
no lag motion occurs. 

Among the parameters presented in ([2]), the stability is usually studied 
with respect to the two parameters which have the greatest influence on it. 
These parameters are <7, which is the distance of the center of gravity (c.g.) 
of the blade from the center of pressure, expressed as a fraction of the chord; 
and the non rotating torsional natural frequency u \. The other parameters, 
described in table [TJ are taken as fixed. Notice that the natural frequencies 
v\ and Ai are given as fractions of fi. Table [l] also shows the values used 
for the fixed parameters in the numerical analysis section. These values are 


taken from the work of Faragher (1996). 


Lead or lag 



Figure 1: Typical hinge arrangement and degrees of freedom of a helicopter 
rotor blade. (Source: Wikimedia Commons) 
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Using two sticks in the cockpit the pilot can exert a control moment 
which directly affects the pitch motion of the rotor. The basic effect of a 
pitch change is a change of the lift force produced by the blade. If all blades 
change the pitch simultaneously, the lift of the whole rotor is changed and 
vertical motion ensues. This is the so called collective pitch control and is 
changed using one stick in the cockpit. If the pitch of the blades changes 
as a function of its angular position, a directional force can be created and 
horizontal flight (forwards, rearwards, sideways) is attained. This is the 
cyclic pitch control and is controlled by a second stick in the cockpit. The 
uncontrolled model ([!]) is then expanded using an input signal which 

corresponds to the motion of the lower end of an actuation rod. Following 
again Faragher ([1996]), the controlled model with the parametric values under 
consideration becomes 


M 


m 


+ c 


m 

m 


+ K 




= Bu(ip) = 


3777/vf 

0 


( 3 ) 


Equations ([T]) and ([3]) are models for hovering or vertical flight only. In 
this situation, the matrices in ([2] are constant, and we deal with linear, 
time-invariant (LTI) systems. When forward (or sideways) motion of the 
helicopter is considered, the aerodynamic forces and moments become time 
variant. More specifically, in these cases the matrices in ([2]) are periodic, 
with a time period T = 2irQ or simply T = 2ir when the azimuth angle i/i is 
used as independent variable. This case is certainly more complicated and 
is left aside for future studies. 


Table 1: Parameters of the Dynamic Model 


Par am. 

Description 

Value 

r 9 

Span wise location of blade’s eg 

4.1 nr 

Ch 

Chord 

0.527 m 

1 

Lock Number 

6.95 

Ai 

First flap natural frequency 

1.1 

a 

Pos. of the eg from the center of pressure 

[0,0.08] 


Torsional natural frequency 

[0,\/8] 
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3 Stability Analysis of the Uncontrolled System 


With the state vector x (?/>) = 


9 pep 


.1 T 


state space as the fourth order system 
= A u x(fl>) = 


, equation (JTJ can be represented in 

( 4 ) 


O 22 h 2 

—M~ l K —M~ 1 C 


with O 22 being a 2-by-2 matrix with all elements equal to zero and I 22 the 
identity matrix of the same dimension. 

As the parameters a and vf change while keeping the other constant, 
two types of instability can appear. When A u has a purely real unstable 
root, the unstable mode is known as pitch divergence and it is evidenced 
by a steady increase in the pitch and flap angles. When the unstable roots 
are complex conjugate, the oscillatory instability known as pitch-flap flutter 
appears. 

The boundary at which the transition from stability to pitch divergence 
occurs is given by the presence of a characteristic root at the origin of the 
complex plane. A sufficient condition for A u to have a root at the origin is 
that the determinant of the stiffness matrix is equal to zero, i.e., \K\ = 0. 
From here we obtain 

2 _ 3-fr g a 
1 2c h A? 

To establish the boundary between stability and flutter, it is necessary to 
find a conjugated pair of purely imaginary characteristic roots at a frequency 
lO f. This is achieved by solving the characteristic equation \jufI 44 — A u \ = 0, 
leading to the following two equations in terms of the crossing frequency 

v\ = 2c o 2 j — (6a) 

Ch {fl 2 ^ 2 + 64cu| — 128wjAf + 64Af) 

96r 9 7 (pv'j - 1 J 

The two stability boundaries are presented in Fig. [2] The red, dashed line 
depicts the boundary at which an unstable pole introduces pitch divergence, 
whereas the blue, solid line indicates a change in stability to the pitch-flap 
flutter. These two lines divide the parametric space in four zones. The 
region without shading corresponds to stable parametric selections. The 
light shaded zone at the bottom of the plot (zone I) indicates the presence 
of a pitch divergence mode. The darker shaded region at the top right (zone 
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Figure 2: Pitch-flap stability boundaries for an uncontrolled system. The 
non shaded region represents stable parametric combinations. The shaded 
regions indicate pitch divergence (I), pitch-flap flutter (II) and both (III) 
unstable behaviors. 


II) indicates pitch-flap flutter behavior, and the darkest region to the right 
(zone III) corresponds to the worst possible conditions, when both unstable 
modes are present. 


4 Stability Analysis of the Controlled System 


Faragher ([1996]) and |Krodkiewski and Faragher |2000) present a method to 
increase the size of the stable region in Fig. [2] by means of delayed feedback. 
Based on the work of Pyragas (1992), the authors of these works propose 
the control law 


u(ijj) = a (/? (t/0 - P (ip ~ t)) + b (ft {ij}) - P (^ - r 


(7) 


with r = 27T. The idea behind Q is to use as error signal the difference 


between the current state and its value one period earlier. Faragher (1996) 


and Krodkiewski and Faragher (]2000) claim that making the delay exactly 
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equal to the period of the motion greatly increases the region of asymptotic 
stability of the uncontrolled system. On the contrary, this particular selec¬ 
tion of the delay value actually introduces instability into the system. In the 
following paragraphs, we analyze the effect of the delay on the stability of the 
system and show under which conditions the delayed feedback is beneficial 
or not. 

Under ([7j the controlled system Q is expressed as 

x (ip) = Ax (ip) + A d x (ip — r) (8) 


with 


A = 

E = 


O 22 

-M-'K + E 

0 3777 via ' 

0 0 


I 22 

—M~ l C + F 

0 3777vfb 

0 0 


Ad 


F = 


O 22 O 22 

—E -F ’ 


The stability analysis of a Linear Time Invariant, Time Delay System 
(LTI-TDS) like ([8j can be done following the two propositions of the Cluster 
Treatment of Characteristic Roots presented by Olgac and Sipahi (2002). 
These two propositions state the following: 


• Although an LTI-TDS has infinitely many characteristic roots, which 
change as the delay does, there is only a finite number of possible purely 
imaginary characteristic roots, which appear at periodically spaced 
values of the delay. That is, if s = ju> c is a characteristic root of 
the system for t = r c , the same imaginary root appears when t = 
t c + 2kir/ou c , k = 1, 2, 3 .... 

• When the delay increases an imaginary characteristic root can move 
to the left or to the right of the complex plane, stabilizing or desta¬ 
bilizing the system. This transition is invariant for roots created by 
the periodicity. That is, if by increasing the delay from r c to r c + e, 
0 < e « 1 , the root at s = juo c moves to the left, it does the same 
when the delay increases from r = r c + 2 tt/oj c to r = r c + 2 i t/uj c + e. 

According to these ideas, to determine the values of the delay for which 
the system is stable we must find the purely imaginary characteristic roots 
of system ([ 8 |. For this, we form the characteristic equation 

det ( 5/44 — A — A d e~ rs ) = 0. (9) 
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Since Ad has rank 1, this equation has the form 

P(s) + Q(s)e~ TS = 0. 


(10) 


By replacing s = joo c in (|10[) , and reorganizing we obtain 


P{Wc) 
Q (^c) 


— e~ jUcT 


( 11 ) 


After equating the magnitudes of both sides of © we arrive to an expression 
in terms of uj c 

\P (u c )\ 2 - \Q (0Jc)\ 2 = o. ( 12 ) 

For the particular case at hand, equation ( |12[ ) is an eighth degree polynomial 
in ui c , with nonzero coefficients only for the even exponents of u c . Introducing 
a new variable rj = we get an expression of the form 


V 4 + C 3V 3 + C2 rj 2 + C\T] + c 0 = 0. 


(13) 


The solutions of (13) with rj > 0 correspond to valid crossing frequencies 
oj c = A^/rj for the LTI-TDS system ([8|. This equation shows that there can 
be up to 4 crossing points at the imaginary axis. To find the delay at which 
such crossing frequencies occur, we equate the phase angles of both members 
of © and solve for r c 


r = — (ZQ (lo c ) - ZP (u c ) + 7r). 
UJ C 


(14) 


The explicit versions of (13) and (14) are very long, and thus left outside 
of this document. 

To find whether an imaginary root crosses to the left or to the right of the 
complex plane, the root tendency for each cu c and one of its corresponding r 


is found using its definition (Olgac and Sipahi, 2002) 


RT = signurn ( 5fi 


ds 

dr 


(15) 


The root sensitivity ds/dr is found from Q using implicit differentiation. 
When RT = 1, the root crosses to the right half of the complex plane, which 
introduces instability. If RT = — 1 the root moves to the stable left half of 
the complex plane. 


With the help of (13), (14) and (15), we can determine for which values 
of the delay the system is stable given a certain parametric selection. The 
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procedure is simple. After finding the crossing frequencies and the delays 
that generate them we create a table of stability transition points. Starting 
with the number of unstable characteristic roots of the non delayed system 
(which is equivalent to the number of eigenvalues of A with positive real 
part), we increase r until a crossing delay is reached. If RT = 1 for that 
point, we increase the number of unstable roots by 2, or decrease it if RT = 
— 1. Those intervals of r for which the number of unstable roots is zero are 
the delays that guarantee stable operation. 


4.1 Stabilization of Pitch Divergence 


Faragher (1996) states that the control action ([Tj) cannot stabilize a system 
that presents pitch divergence. They arrived to this conclusion by means of 
a parametric search. Here we provide an analytic proof of it. 

When the helicopter rotor presents pitch divergence, the non delayed 
case has a real positive characteristic root. For this root to become stable, 
a stabilizing transition at the origin of the complex plane must take place. 
That is, we need a crossing frequency u = 0 for which RT = — 1. For (13) 


to have a root at uj = 0, the independent term Co must be zero. Considering 
the explicit form of this coefficient, the condition for this to happen is 


Co = (Af — 944ja) uf + ( 3^2 _ 14167 a) v\ + 

Cfi ^ 


97 2 r 2 g o 2 


= 0. (16) 


This equation shows that for a, isf > 0 it is impossible to set the control pa¬ 
rameters to a value such that Co = 0. Therefore no crossing of the imaginary 
axis occurs at the origin (in either direction) and the transition from pitch 
divergence to stability cannot take place. 


4.2 Stabilization of Pitch-Flap Flutter 

In order to show that the proposed control law 0 is able to eliminate the 
pitch flap instability in a wide range of values of the operating parameters, 
explores the complex planc 0 s = % + jcn searching for lo¬ 
cations where the loci of zero points of the real and imaginary parts of the 
characteristic equation ([9]) intersect. Their variation of parameters follows 
this logic: after selecting the values of a and 17 , the control gains a and b 
are varied, searching for the optimal location of the characteristic roots, i.e., 
the location where the rightmost roots have the smallest possible real part. 

x We use x — K(s) to avoid confusion with the parameter a. 


Faragher (1996 
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Figure 3: Root locus of zero points of the real (blue, solid) and imaginary 
(red, dashed) parts of the characteristic equation ([9]) for a = 0.08, uf = 10.8, 
a = 6.75 x 10 , b = 0.6 x 10 -4 , t = r. The intersections indicate the 

location of characteristic roots of the system. 


This method has a big drawback: since (J9]) is a transcendent equation 
it has infinitely many roots and any parametric exploration would never be 
able to find all of them to determine which one is the rightmost. We used the 
same searching method they applied and did not obtain the same root loci 


and characteristic roots presented by Krodkiewski and Faragher (2000) in 


figure 5 and by Faragher (1996) in figure 6.19. For the parametric selections 


stated in those workqjand considering the same search region, we obtained 
the plot in Fig. 3(a)| The presence of three unstable characteristic roots 
is clear. Furthermore, one of these roots is real, implying that the control 
action is introducing pitch divergence into a parametric selection that did 


not exhibit this behavior. The roots in Fig. 3(a) are, however, not rightmost. 
If the search region is increased further to the right, another real root shows 
up, as it is shown in Fig. |3(b)| This highlights how parametric searches of 
this kind may miss very important information. 

In order to verify the results of Fig. |3(b)[ we used the QPmR function 
(Vyhlfdal and Zftek 2009) to numerically find the roots in the same region 


of the complex plane. The results coincide: QPmR finds the roots at 0.9927, 
0.0595, -8 x 10“ 4 + jl.0623, -1.15 x 10” 2 + j2.0439, 4.10 x 10“ 2 + j2.8268 
and —0.196 + j3.5967. We also simulated the system using the same initial 
condition stated by Krodkiewski and Faragher (2000): [9 /3] T = [0 0.01] T . 


2 Please note that in this paper we are not using the 10 4 scaling factor in the control 


gains introduced by Faragher (19961 or Krodkiewski and Faragher (2000). 
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(a) T = 2n 



(b) T = 0.2296 


Figure 4: Simulation results for a system with a = 0.08, uf = 10.8, a = 
6.75 x 10~ 4 , b = 0.6 x 10~ 4 with two different delays. 


Table 2: Crossing Frequencies and Delays for a = 0.08, v‘\ = 10.8, a = 
6.75 x 10 -4 , b = 0.6 x 10~ 4 _ 


^ Cl 

RT 

7"cz0 

2.1949 

-1 

0.0852 

1.0525 

1 

0.3580 

3.0268 

1 

1.519 


The simulation results of Fig. 4(a) show a time constant which is consistent 
with the presence of a real unstable root located near to s = 1. 

For the selected control parameters in the previous example, which are de¬ 
clared as the optimal, the delay value selected by[Krodkiewski and Faragher 


(2000) actually creates instability. By applying the CTCR methodology de¬ 
scribed at the beginning of the section, we can find which values of the delay 
are actually able to stabilize the system. For this particular parametric se¬ 
lection, the solution of (13) gives three crossing frequencies, which are listed 
in Table [2] together with their corresponding delays and root tendencies. 

The smallest crossing delay has RT = — 1. This implies that the system, 
which presents pitch-flap flutter instability for r = 0, recovers its stability 
when the delay is larger than r = 0.0852. However, it loses its stability 
for delays larger than r = 0.3560. The optimal delay value is then in the 
interval [0.0852 0.3560]. By varying the delay whilst keeping the control 
gains fixed, which is a numerical search approach similar to that used by 
Krodkiewski and Faragher (2000|, we found the optimal delay to be r = 


0.2296, as it is shown in Fig. [5] The rightmost characteristic roots for this 


11 




































| -Real Part of the CE-Imaginary Part of the CE | 


Figure 5: Root locus of the real (blue, solid) and imaginary (red, dashed) 
parts of the characteristic equation ([9]) for a = 0.08, uf = 10.8, a = 6.75 x 
10 -4 , b = 0.6 x 10 -4 , t = 0.2296. The intersections indicate the location of 
characteristic roots of the system. 

case are —0.4368 + jl. 2018 and —0.4368 + jl. 9596. Figure [4(b)| presents 
the simulation results for this case. The search of the optimal delay value is 
made possible by the knowledge of the exact stability interval with respect 
to the delay. Knowing the upper and lower bounds, we can search in a finite 
space where we are sure to find the rightmost root with smallest possible 
real part. 

The optimization just presented assumed a fixed set of control gains and 
changed only the delay. By varying the three parameters it may be possible 
to move the rightmost root further to the left to increase the stability margin 
of the system. Figure [6j for example, presents the real part of the rightmost 
root of the system as a function of the control parameters a and b for a fixed 
delay value of t = 0.2296. Here we observe that the control gain a has a 
stronger effect than b on the location of the rightmost root, and that the 
largest stability margin for this delay is obtained with control gains slightly 
larger than those selected. 

A better estimation of the optimal control gains for this delay value 
is made possible by observing Fig. [7J the contour plot of the surface in 
Fig. m For example, the control gains a = 7 x 10 4 and b = 1.03 X 10 4 , 
marked by a black square, locate the rightmost root of the system at s = 
—0.4475 + j’1.7690, giving a settling time a bit smaller than the one obtained 
with the current selection of control gains, marked with a blue dot in Fig. [7] 
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Figure 6: Real part of the rightmost root of the system as a function of a 
and b for a fixed delay value of r = 0.2296. 


The knowledge of the stability regions in the domain of the delays allows 
the user to search for the best delay given a certain control gain selection. 
Since we are dealing with retarded systems, the stability region does not 
change dramatically as the parameters change and therefore a gradient de¬ 
scent algorithm can be implemented. With such systematic optimization 
routine, the best possible parametric (a, b , r) combination for a given set of 
(<r, v\) parameters can be found using this approach. 

5 Closing Remarks 

In this document we studied the problem of stabilization of pitch-flap mo¬ 
tion of helicopter rotors using delayed feedback. We revisited a previously 
proposed control law and by using newer tools for the analysis of linear 
time-invariant, time delay systems we showed how the delays which make 
the system stable can be found. With his knowledge it is possible to find the 
largest possible stability margin of the system and we suggested a method¬ 
ology to do so. 

The Maple and Matlab codes used to create the plots in this document 
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Figure 7: Contour plot of the surface in figure [6] 


can be downloaded from https://db.tt/MHvoLtdM. 

This document considered only the case of hovering flight, for which the 
system matrices are constant. In the case of forward flight, the dynamic 
model of the system becomes time periodic, and thus different tools need to 
be used for its analysis. This topic is a matter of further research. 
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